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Abstract 

We reconsider the problem of calculating arbitrary negative integer moments of the (regu- 
larized) characteristic polynomial for N xN random matrices taken from the Gaussian Unitary 
Ensemble (GUE) . A very compact and convenient integral representation is found via the use 
of a matrix integral close to that considered by Ingham and Siegel. We find the asymptotic 
expression for the discussed moments in the limit of large A'^. The latter is of interest be- 
cause of a conjectured relation to properties of the Riemann function zeroes. Our method 
reveals a striking similarity between the structure of the negative and positive integer mo- 
ments which is usually obscured by the use of the Hubbard-Stratonovich transformation. This 
sheds a new light on "bosonic" versus "fermionic" replica trick and has some implications for 
the supersymmetry method. We briefiy discuss the case of the chiral GUE model from that 
perspective. 



1 Introduction 

Recently there was an outburst of research activity related to investigating the moments and 

correlation functions of characteristic polynomials Zn{^) — det (^(iIn — for random N x N 

matrices H of various types [Q, |, |, |, |, |, 0, |, |, |l^, 

There are several, not completely independent, sources of motivation behind studying character- 
istic polynomials. First is the intriguing conjecture relating limiting distribution of the non-trivial 
zeroes Sk = \ + itk of the Riemann zeta function C(s), on the scale of their mean spacing, to that of 
(unimodular) eigenvalues of large random unitary matrices . This implies that locally-determined 
statistical properties of C(s)j high up the critical line Res — 1/2, might be modelled by the corre- 
sponding properties of 2'(/i), averaged over the so-called Circular Unitary Ensembles (CUE), i.e. 
with respect to the normalized Haar measure of the group U{N) oi N x N unitary matrices. Such 
a line of thought and underlying evidences in favour of the conjecture are explained in detail in 
the papers by Keating and collaborators, see [0, |[ ||]- In particular, in the authors managed to 
evaluate arbitary moments of \Zn {fJ- — e'^) | explicitly: 

/\7 f •i0M2"\ -rr r(i)r(j + 2n) MN.2n j! . 

where T{z) is the Euler gamma-function. The moments as above were derived for Re n > —1/2 
but can be analytically continued to the whole complex n— plane. The limiting value is presented 
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for the integer positive n. This should be compared with the conjecture ||l|, 

-I rft|C(l/2 + .i)r^a„jn(-^(log-r^ (2) 

for the values of the positive integer moments of the Riemann fmiction as T — )■ oo. Here a„ is 
the number specific for C-function[Q|, but the rest shows universal features common to both the 
random matrix calculations and function. The parameter 57 log (^) in the above equation 
plays the role of the inverse spacing between the (— function zeroes at a height T and should be 
identified with N/2tt of the unitary random matrix calculations 0] . 

It is important to have in mind a high degree of universality of the obtained results, as discussed 
in the work by Brezin and Hikami[Q . By universality one usually means insensitivity of the spectral 
characteristics to details of distributions of matrix entries. In particular, the limiting value Eq.(l) 
of the positive integer moments of the characteristic polynomials for unitary matrices is shared, 
after an appropriate normalisation, by a broad class of Hermitian random matrices, whose most 
prominent representative is the Gaussian Unitary Ensemble (GUE). 

Other quantities like the distribution of the logarithm of the characteristic polynomial, its 
derivative, etc. enjoyed thorough investigations as well|^, ||, |[ ||. The results were also extended 
to the ensemble of unitary symmetric matrices (COE- Circular Orthogonal Ensemble), which are 
related to statistics of zeroes of the so-called L— functions |^|. An updated summary of open 
questions on relations between the properties of Riemann zeta function and random matrices can 
be found at web-page of the American Institute of Mathematics pO| . 

What concerns negative moments of the characteristic polynomials, an additional interest in 
calculating them arose because of a conjectured behaviour of the negative moments of the (regu- 
larized) Riemann zeta function: 
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put forward in Q for 1 < 5 < logT. 

The formula Eq.(|l|) shows divergency at negative integers n and thus provides one with no 
explicit answer. Such a divergency is a natural consequence of necessity to regularize characteristic 
polynomials by adding a small imaginary part to the spectral parameter /i to avoid singularities 
due to eigenvalues. When such an imaginary part is comparable with the separation between the 
neighbouring eigenvalues one again might expect universality of the corresponding expressions. 

The Section 8 of the work by Brezin and Hikami [Q discusses a possible way of calculating 
the negative moments of the characteristic polynomials themselves. However, in contrast to the 
moments of the absolute values those are not divergent and, when taken alone, are insufficient for 
the sake of comparison with Eq.(||). 

Original goal of the present paper was to reconsider the problem of calculating both the negative 
integer moments of the characteristic polynomials and those of their absolute value. We succeeded 
in the analysis of correlation functions of the (regularized) characteristic polynomials in the limit 
of large N and obtained: 
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where the regularizations Im/ii > 0,lnifi2 > as well as the spectral difference u = Re(/ii — 112) 
were considered to be of the order of mean eigenvalue spacing = [Np{p))~^], with p{p) being 
the mean eigenvalue density at /x = ^Re{^i + fi2). 

Such an expression complements that given in Eq.(|^) and is expected to be universal and 
applicable to the Riemann zeta-function. Indeed, taking into account the nonuniformity of the 
spectral density for GUE the correspondense between the parameters should be as follows: A^^ — 
Np{fi) ^ log (T/2tt). We see that the random matrix result Eq.(^ and the conjectured Riemann 
(■-function behaviour Eq.(H) agree in the overall parametric dependence. 

Another motivation for such a calculation comes from some questions that arose in applications 
of the random matrix theory to chaotic and disordered quantum systems which we shortly discuss 
below. 

As is well known, eigenvalues of large random matrices pl[ ^ played a prominent role in 
the development of the field of quantum chaos, see e.g.||2^, and in revealing its connections to 
mesocopic systems as well as to some aspects of Quantum Chromodynamics [p5| . Results 
on moments and correlation functions of the characteristic polynomials of large random matrices 
related to various aspects of quantum chaotic systems can be found in @ 0, |, §, see also 
p^ , p7[ | for related studies. Characteristic polynomials for chiral random matrix ensembles are 
used as a model partition function for the phenomenon of chiral symmetry breaking and as such 
enjoyed thorough considerations, see and references therein. 

Intimately connected with the field of quantum chaos is the domain of mesoscopic disordered 
systems. The paradigmatic example is a single non-relativistic quantum particle moving at zero 
temperature in a static random potential. The system Hamiltonian is, in essence, equivalent 
to a matrix with random entries. Moreover, in some limiting case such a matrix belongs to 
the "domain of universality" of the classical random matrix theory. This fact is of paramount 
importance and follows from the seminal Efetov's work, see the book where the notion of 
the supermatrix (graded) non-linear cr— model was introduced for the first time. The latter tool 
alternative to other techniques in the theory of random matrices expresses expectation values of the 
(products of) resolvents of random operators in terms of integrals over graded matrices containing 
both commuting and anticommuting entries. The method proved to be capable of dealing with 
quantities less accessible by other methods and turned out to be indispensable in establishing links 



between the theory of random matrices and quantum chaotic/mesoscopic systems, see|26, 2^ and 
references therein. 

An alternative technique which enjoyed many applications in theoretical physics of disordered 
systems of interacting particles is the (in)famous "replica trick" . Suppose one likes to calculate the 
ensemble average (log Z) of a logarithm of some quantity Z. The replica trick exploits the relation: 
logZ = lim,i^o -(•^" — 1) and attempts to extract the averaged logarithm from the behaviour of 
the moments {Z"), with n being either positive or negative integer. It is clear that in general the 
limiting procedure suffers from non- uniqueness of the analytical continuation and " mathematicians 
will throw up their hands in horror or despair, while physicists are much intrigued" [p9| . Random 
matrices provide an important testing ground for the replica calculations, with the absolute value 
\Zisf{p,)\ of the characteristic polynomial playing the role of Z. The advantage here is that one has 
a better control on results obtained by the ill-defined recipe comparing them against those known 
from independent calculations. 

In particular, the early paper by Verbaarschot and Zirnbauer devoted to the relation 
between the replica and supermatrix methods revealed inherent problems in the former absent in 
the latter. They found that the natural analytic continuation n gave two different answers for 
the "fermionic" (positive moments) and "bosonic" (negative moments") versions of the replica trick. 



3 



neither of them coinciding with the known result. In contrast, the latter is correctly reproduced 
within the supermatrix approach. 

Very recently the verdict of inadequacy of the fermionic replicas was challenged by Kamenev 
and Mezard||ll| and further elaborated by Yurkevich and Lerner||l^. In particular, Kamenev and 
Mezard discovered a convenient integral representation for the integer positive moments providing 
one with a better control on analytical structure of the expressions. This allowed them to put 
forward an ansatz which yielded in the limit n — > the correct exact (nonperturbative) result 
for the GUE matrices, and the correct asymptotic results for other symmetry classes . A critical 
analysis by Zirnbauer p9| ] demonstrated in a coherent manner that the proposed ansatz was in 
no way a well-behaved analytical continuation. Even so, such a critique did not devaluate the 
recipe itself but rather restricted its domain of applicability to perturbative calculations and called 
for further investigations. And indeed, the amended fermionic replica trick immediately found 
applications in the theory of disordered electronic systems with interactions ^ when it was 
among very few tools actually available. Let us also mention a recent development in the framework 
of the Calogero- Sutherland model inspired by closely related ideas 

The discussed new insights in the nature of the fermionic replica left, however, unclear if one 
could come forward with a meaningful amendment for their bosonic counterpart within the context 
of nonlinear cr-model ideas (see, however, ||30|] for the replica limit in the context of orthogonal 
polynomials). 

An additional motivation for the present paper was to try to bridge the gap between the cases 
of the positive and negative n. Our attempt succeded in discovering an integral representation for 
the negative integer moments which is strikingly close to that obtained by Kamenev and Mezard 
for the positive ones. 

Technically, analyticity properties inherent in the negative moments of the absolute value of 
characteristic polynomials is known to result in the non-compact (" hyperbolic" ) nature of the in- 
tegration manifold for the bosonic nonlinear a— model discovered by Schafer and Wegnerjs^. In 
standard considerations such a manifold enters via the so-called Hubbard-Stratonovich transfor- 
mation (see the Appendix D for more details). It came as quite a surprise to the present author 
that the Hubbard-Stratonovich transformation turned out to be not only unnecessary, but played, 
in fact, a misleading role hiding the simple structure of the negative moments. To reveal that 
structure one should introduce an alternative route via use of the matrix integral close to one 
considered by Ingham[|6| and Siegel[Q many years ago. 

As to the replica limit, the fact of close similarity between our integral representation and those 
in ^ makes it apparent that very the same KMYL recipe "works" for the bosonic version in 
the same way as for its fermionic counterpart. This should not be considered as contradicting the 
Zirnbauer's argumentation since both versions of the replica trick are somewhat deficient, in the 
strict mathematical sense. The result obtained just indicates that accepting one of them we have 
little reasons for discarding the other. 

Clearly, our way of dealing with negative moments suggests certain revision of the supermatrix 
method whose underlying technical idea is a simultaneous uniform treatment of both types of the 
moments (positive and negative). In fact, we show that after the disorder average is performed 
treating "fermionic" and "bosonic" sectors differently can be of some advantage. 

The structure of the paper is as follows. In the section II we expose our method on the simplest 
example of negative integer moments of the characteristic polynomials and analyse the obtained 
expressions in the limit N ^ oo. Then in the section III we proceed through the calculation for 
the negative moments of the absolute value of the polynomial (in fact, a correlation function). In 
the section IV we comment on the replica trick and illustrate our statements by addressing briefly 
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the case considered in |33 - the chiral GUE model - from that perspective. Finally, in the section 
V we present the simplest nontrivial example of extention of our method to the general type of the 
correlation (generating) function containing simultaneously both positive and negative moments of 
the characteristic polynomials of GUE/chiral GUE matrices. The open questions are summarized 
in the Conclusion. Technical details are presented in the appendices. 



2 Negative Moments of the Characteristic Polynomial 

Let H he N X N random Hermitian matrix with characterized by the standard (GUE) joint 
probability density: 

V{H) = Cr^exp-jTTH', C^r = (27r)-^^7V^V2 (5) 

with respect to the measure dH — YiiLi dHu Yli<j dHijdH* ij . Here we use * to denote complex 
conjugation and denote: dzdz* = 2c?Rezc?Imz. 

Regularizing the characteristic polynomial Zi^{fi) — det (^(iIn — by considering the spectral 

parameter /i such that Im/i > one represents negative integer powers of the determinant as the 
Gaussian integral: 

„ n f . n .71 ^ 

[^A^(A^)""]-7i-^^/ n^'S.expa^^StS.-^^StTJS. 

^ ' fc=i I fc=i fc=i J 

where for k — l,2,...,n we introduced complex A^— dimensional vectors = (sfc.i, Sfc.w)"^ 
so that d^Sk = Y[iLidsk,ids'l ^ and T, f stand for the transposition and Hermitian conjugation, 
respectively. 

Denoting by (...) the expectation value with respect to the distribution Eq.(|^) we are interested 
in calculating the negative integer moments of the two types: 

^SnlA^l) = ([^iv(/il)]-") (7) 

as well as 

/cj^;„(A*i>M2) = {[ZN{f^i)ZNil^*2)r") (8) 

assuming the regularization Im(/ii) — Im(/i2) > 0. In particular, when Re/ii = Re/i2, the latter 
quantity amounts to the negative moment of the absolute value of the characteristic polynomial. 

Let us start our consideration with the simplest of the two. Performing the ensemble averaging 
in the standard way one finds for the moments of the first type: 

- 7i-wv / n d'Sk exp 1;., Y: sIs^- - ^ E (sis.) (S^S,) (9) 

^ ^ k=l k=l k,l=l J 

Further introducing a n x n Hermitian matrix Q with the matrix elements Qki — S|S; the 
integrand is conveniently rewritten as: 

expj^/iiTrQ-^TrQ^I 
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The standard trick suggested to deal with the apparent problem of the non-Gaussian integral 
above is to employ the famous Hubbard-Stratonovich transformation amounting to: 



exp 



-^TrQ^I ^ f dQexp |-f TrQ^ _ IxrQQ 



(10) 



thus trading the integration over n x n Hermitian matrices Q for a possibility to perform the 
Gaussian integration over the vectors Sfc. Then the resulting matrix integral is amenable to the 
saddle-point treatment in the limit N oo. 

However one may notice a possibility of an alternative route. Its starting point is similar to 
the method employed in ^8, 39 1 where it was suggested to rewrite the integral Eq.(^ introducing 
the matrix (5— distribution as the product of (5-distributions of all relevant matrix elements. Then, 
obviously, 

- f dQe-^TrQ^/„(Q) 



JC 



where 



IniQ) 



k=l 



OC 



S, Si, 



k<l 



SIS; 



(11) 



(12) 



and the (5-distribution for complex variables is understood as the product of the (5-distributions 
for their real and imaginary parts. From now on we do not take care explicitly of multiplicative 
constants in front of the integrals. We will show how to restore the constants on a later stage using 
the normalisation condition. 

To evaluate the last expression we employ the Fourier integral representation for each of the 
delta-functions involved and combine the Fourier variables into a single n x n Hermitian matrix 
F. This allows us to proceed as follows: 



IniQ) OC / nd^g^g-MiELi^IS'- / dF exp l'-Tr(^FQ)-'-Y,Fik {sis 

•' k=l •' y. kl 



(13) 



dPe^^'^'y^^) det(F-/iilr 



-N 



Up to this point our consideration was, in fact, parallel to that employed in js^, We 
however suggest to go one step further by noticing that the last matrix integral is quite close to 
the distinguished one considered originally by lngham|36 and Siegel|37j: |^ 



di^e-Tr(FQ) 



F>0 



detF 



(27r)- 



-]5!(p+l)!...(p + n- l)!detQ-(''+") (14) 



where both F and ReQ are positive definite Hermitian of the size n and the formula is valid for 
p > 0. The Ingham-Siegel integral can be viewed as a direct generalisation of the Euler gamma- 
function integral: T{j) + l)g~'^P+^) = /^^q dffPe~^'^ to the Hermitian matrix argument and paved 
a way to the theory of special functions of matrix arguments which is nowadays an active field of 
research in mathematics and statistics, see e.g pO||. 



^In fact, Ingham and Siegel considered the set of real symmetric matrices F rather than their Hermitian coun- 
terparts and found the result: (tt) J Ilfc— i ^ ^^2^) '^'^^ "'^^^'''^ ~^ ■ However, their method is equally 
applicable to both cases. 
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It is an easy matter to adopt their method to calculating our integral |^ which is a matrix- 



argument generalisation of the formula: df 



^i^q)^-^e'f^' for g > and zero oth- 



erwise, provided Im/i > 0. Performing the calculation (Appendix A) we find for N > n: 



In,N{Q > 0) = / dPe^^^^^'^) [det (f - mi1„ 



-N 



= CAr^„detQ^-"e5^iTrQ 



(15) 



with C 



Na 



(2^)' 



and In,N{Q) — whenever at least one of the eigenvalues of Q is 



negative (we recall our choice Im/ii > 0). 

As a result we arrive (after rescaling the integration variable: Q 2NQ) to the following 
integral representation for the negative integer moments of the characteristic polynomial in terms 
of the integral over the matrices Q: 



(1) 



c 



(1) 



f dQe-w[-'wTrQ+iTrQ^] qN-u 
Jq>o 



(16) 



provided N > n. 

The overall constant c'^\ can be restored by noticing that for Re/zi oo the moments tend 
asymptotically to On the other hand, it is easy to understand that such a limit is equivalent 

to discarding the quadratic in Q term in the exponent of Eq. ( [T^ ) . The resulting integral is precisely 
the Ingham-Siegel one, Eq.(n4r), and comparison yields the required constant: 



C 



(1) 

N,n 



{-iNy^'\2Tr) 



nN-l ., 



As the last step of the procedure we choose eigenvalues gi, <?„ and the corresponding eigen- 
vectors of (positive definite) Hermitian matrix Q as new integration variables. This corresponds to 
the change of the volume element as: dQ ~ GnA^jg} Yl7=i dqidii{Un) where the factor A'^jg} — 



ni<j(9i ~ squared Vandermonde determinant, G„ = (Stt)" ^ ypr — ^ 

stands for the normalized invariant measure on the unitary group U(n). The integrand is obvi- 
ously U (n) invariant and we obtain: 

det (/ii 1 AT — H) 



and dfj,{Un) 



JV.i 



qi>0 



]l{dq,qr)A'{q}exp-NJ2AiQ^) (17) 



where 



C 



(1) 

N,n 



Nr, 



and A{q) 



z/iig — In g 



(18) 



The last integral representation is our main result for the negative moments of the first type: 
^Ar^n(Mi)' valid for arbitary N > n. One can further play with the formulae for finite N and n, 
expressing, for example, the negative moments as n x n determinants: 



det (^ilAf — H) 



cx det [^jk] 



(19) 



d<?'?^-'V(.i)(g)4^)(g)e^[^'^^''-i«^ 



^We suggest to call such an integral "the Ingham-Siegel integral of second type". 
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where '!Tj^\q),TTj^\q) are any monic polynomials of degree j in the variable q, compare with the 
case of positive moments in . 

In practice, however, we are mostly interested in the limit of large matrix sizes where one 
expects the results to show universality as was discussed in much detail in the Introduction. To 
extract the leading asymptotics as — > oo when keeping moment order n fixed one should employ 
the saddle-point method and find the saddle points of A{qi). 

Before doing this we observe that the structure of the derived expressions show striking simi- 
larity to those obtained for the positive moments of the characteristic polynomials found in [pl| , 
see also |, |l§: 

det(Miljv-i?) ) = C'^'le^'^M l[dq,A^{q}e^p-NJ2MQ^) (20) 



= 1 



where 



and the expression for A{q) is the same as in Eq.(|]i 

The only essential difference between the two representations (apart from that in the multi- 
plicative constants and a slight change of the power of the determinant: N — n rather than just 
N, which is anyway irrelevant for large N) is the range of integration. For the positive moments 
one integrates over the whole real axis — oo < < oo whereas it is over the positive semiaxis 
< qi < oo for the negative moments. 

Thus, we need to consider the saddle points of A±{q). It is convenient for further reference to 
define fii = fi + ^ + iS, with /i, oj, S -real, and consider Nuj, N5 to be fixed when N oo. Then 
one can replace fii with /i in the saddle-point calculations. The saddle points are obviously given 
by equations: 

q^-i|^--=0 (21) 
where i — 1,2, n. Each of these equations has two solutions: 



,± = ^A^±V4^ (22) 

We would like to choose the spectral parameter fi to satisfy \fi\ < 2 in accordance with the idea 
of considering the bulk of the spectrum for GUE matrices of large size. Then only for g"*" the real 
parts are positive and the corresponding saddle points contribute to the integral over the positive 
semiaxis: q > 0. Consequently, among 2" possible sets of saddle points (^qf, q^) only the choice 

g+=diag(g+,...,g+) (23) 

should be considered as relevant. This feature constitutes a considerable difference from the case 
of positive moments where all 2" saddle-points yield, in principle, non-trivial contributions, albeit 
of different order of magnitude in powers of the small parameter N~^. For example, for n = 2K 
the leading order contribution in the later case comes from the choice of half of saddle-points to be 

q^, the rest being q^ , with the combinatorial factor ( j counting the number of such sctsQ. 
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Presence of the Vandermonde determinants makes the integrand vanish at the saddle-point sets 
of the exponent and thus care should be taken when calculating the saddle point contribution to 
the integral. This part of the procedure uses explicitly the so-called Selberg integral: 



-°° k=l ki<k2 3 = 1 



(24) 



for i > 0, see the paper by Kamenev and Mezard for more details. General points of their 
analysis are applicable for our case without any modification. 

Expanding around the relevant saddle-points: qk — (1^ S.k and performing the required cal- 
culations we find in a straightforward way the asymptotic expressions for the negative moments: 



(-i)^"iV^" 
I iujNn 



-4- c^^y 



t/2 



■prAf-1 ., 



(4-^2)-^ 



(25) 



X exp ■ 



(zM + V4~J^) 



Nn 



The formula for A^^\j(m2) where = n — ^ — i5 can be obtained from the above expression by 
taking its complex conjugate and changing uj —lu. Taking the product of the two expressions 
we finally find: 



fN-l 



iTrp{fi)uj 



1 



2 



(26) 



where we used the known expression for the (semicircular) mean density of GUE eigenvalues: 

This completes the calculation of the denominator in the formula Eq.(^. To find the corre- 
sponding numerator we proceed to derivation of the analogous expressions for the moments of the 
second type. 



3 Correlation functions for the negative moments of the 
characteristic polynomials. 

To this end, we consider the product of the expression Eq.(|^) with its complex conjugate at a 
different value of the spectral parameter and average it over the GUE probability density. From now 
on we use the index a — 1, 2 to label the N-component vectors S^- stemming from the first/second 
set of the integrals. To write the resulting expression in a compact form it is again convenient to 
introduce 2n x 2n Hermitian matrix Q with the matrix elements Q^; ''^^ = S^^ ^.S^a.;, with k and 
I taking the values 1, n. In terms of such a matrix we have: 



^(2) 



/n f . n .71 ^ '\ 

Yl (fSi.k(fS2,k exp I ^fii S^Si,fc - S^.S2,fc - — Tr (qLQL^ I 

fe=l I fe=l k=l ) 



(27) 
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where L = diag(l„, -1„). 

Again, the standard way is to use a variant of the Hubbard-Stratonovich transformation Eq.(|o|) 
allowing to convert the term quadratic in Q (quartic in S) to that linear in Q (quadratic in S) and 
integrate out the vectors S. However, presence of the matrix L and the requirement of convergency 
of the Gaussian integrals necessitates introducing this time a rather non-trivial domain (the so- 
called "hyperbolic manifold", [Q) for the integration over Q, to make such a "decoupling" well- 
defined. This problem comprehensively discussed e.g. in p8[ |4l| makes the whole procedure 



technically involved. For a good pedagogical introduction sec [p3| , the outline of the procedure is 
presented in the Appendix D of the present paper. 

For the method suggested in the present paper such problem does not arise at all. The 2n x 2n 
matrix Q is a Hermitian positive definite and the whole procedure at this stage does not require 
any modification. Employing the Ingham-Sicgcl integral of second type yields in this case: 

ICPjf^^,H2)^dA / dQe-^HTrA^Q+^Tr(QLQL)]^^^QAr_2„ , m = diag(Mil„, -A^^ln) 

JQ>0 

(28) 

provided N > 2n, with the overall constant 

= /Ar^2^^n/r, ^-n(2n-l). 



[lj=N-2n3- 

Clearly, such a uniform applicability can be considered as a technical advantage. Nevertheless 
hyperbolic structure, in fact, lurks in the expression above and manifests itself at the next stage. 
Namely, equation Eq.(^8|) differs from its analogue Eq. (^6|) in one important aspect: it is now 
of little utility to introduce eigenvalues/eigenvectors of Q as integration variables. Rather, it is 
natural to treat Ql — QL as a new matrix to integrate over. Such (non-Hermitian!) matrices are 
just those forming the mentioned hyperbolic manifold. I find it sensible to discuss their properties 
explicitly in the Appendix B. They satisfy = LQ^L, have all eigenvalues real and can be 
diagonahzed by a (pseudounitary) similarity transformation: Ql = TqT~^, where q = diag((7i, (72), 
and n X n diagonal matrices (71,92 satisfy: qi > , 52 < 0. Pseudounitary matrices T satisfy: 
T'^ LT = L and form the group f7(n, n) ("hyperbolic symmetry"). 

In fact, a more convenient way is rather to block-diagonalize matrices Ql as 

QL=fJP^ ^ Vo-^ , Where To. ^^^'^^ 



P2 / ^ ' U{n) X U{n) 

and Pi^2 are n x n Hermitian, with eigenvalues qi^2, respectively. The integration measure cZQl 
is given in new variables as [^: dQ^ = dPidp2 Hfe^ (^i.fci — 92, ^2)^ d,^(T) where the last factor 
is the invariant measure on the manifold of T— matrices whose explicit expression is presented for 
reference purposes in the Appendix C. 

We therefore arrive to the following expression: 



/cj^L (X / / dPidP2l{M,Pi,P2) (29) 

' Pi>0 J Pi<0 



X n ('?i,..-'?2,.J^detP-~^"det(-P2)"""e--Tr(P^+P, 



fci ,k2 
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where 



I{M, Pi, P2) = / dfi{T) exp jiiVTr 

1 



Tn 



cx [~i{^J.l - Afa)]' 



Pi 



(30) 



The calculation of the above integral is presented in the Appendix C. We see that its value 
depends only on the eigenvalue matrices qi and 92- As a final step we change P2 —P2 and again 
introduce those eigenvalues (and corresponding eigenvectors) of the Hermitian matrices Pi > 
and P2 > as the integration variables. This results in the following expression for the correlation 
function of negative moments of the characteristic polynomial: 



/Cj^j„(^i,/i2) = ( detifiilN~H)detin*2lN-H) 
1 



fi^l/ JO -- JO 

n (,i,., +g2,.je-^i:L.^^(«^^')-^i:L.^^fe.) 

fei ,k2 



n dq,,, q-J- A2{gi} / J] dq,, q'^^ ^^{q,} (31) 



where 



C 



(2) 
N.n 



N 



nf=7v'-jj!(j-n)!] n;Lij! 



and 



-41(9) = - i/^ig - Ing , A2((7) = +iAi29-lng 



(32) 



(33) 



The constant „ given above is most easily checked by considering the limit /ii ^ /i2 ^ 1 in 



IM^nn 

both sides of Eq. (Bll) and using the identity: 

/ Y[dq,q!A'{q}e-^T.l^'^^^p--(-+P)lljl [] 



n p+n—1 



(34) 



valid for p > and Pe/3 > 0. Such a formula is an immediate consequence of Eq.(|l4|) when going 
to eigenvalues of the matrix F as integration variables and considering: Q — /31„. 

Again we see that the structure of the derived expressions is strikingly similar to those obtained 
for the correlation functions of the positive moments of the characteristic polynomials |ll[]: 



det (Aiiljv - H) det (/i^lw - H) 

poo nOC 

^ ^(2)^^-[,? + (,-)^] / Y[dq,,AHq\} / l[dq2,AHq\} 



(35) 



-42(?2,,) 



fci ,k2 
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where 



c 



(1) 



(36) 



the constant C^'^ is defined earUer in Eq.(|20|) and expressions for Ai{q), ^2(9) are the same as for 
the negative moments, Eq.(p2|). 

Now, however, the difference between the domains of integration has more important conse- 
quences. Namely, the negative moments of the absolute value of the characteristic polynomial are 
truly divergent for (/ii — 1JL2) 0, as represented by the factor (/ii — ^2)^" the corespond- 
ing formula. For their positive counterparts such a singularity is fake and is compensated when 
performing the integration along the whole real axis. 

Again, we would like to perform the asymptotic analysis for N — > 00. As discussed in the 
Introduction the most interesting " local" universal regime is to occur when one keeps the difference 
Re(/ii — (12) = ^ and the regularisation 8 so small as to ensure A^max(a;,(5) < cxd in such a limit, 
whereas [i = Re is kept in the range < 2. To shorten our notations we include the 

regularization 5 into w, so that /ii^2 = A* i u jl. Then we can write: 



1=1 



^^+(gi,,) + ^A_(g2,0 



1=1 



4=1 



where the functions A±(ci) are obtained from A\^2{(i) by setting fJ-i = IJ-2 = fJ-- 
The stationary points of A±{q) which are obviously given by the equations: 



91, 



where i = 1, 2, . 



— and (?2 i + 

qi,i ' 12,1 

, n. Each of these two equations has two solutions: 







(37) 



it 



ifj. 



and 



4 



(38) 



but only for q^2 the real parts are positive and the corresponding saddle points contribute to the 
integral over the positive semiaxis: qi^i > or (72,1 > 0. Consequently, among 2^" possible sets of 
stationary point (^qf 



ij ■ 



' 9l,ni 92, 1' 



^2 n) only the choice 
'+ = diag(gj'', ...,q^,q^, ...,qt) 



(39) 



should be considered as relevant. 

Taking care of the Vandermonde determinants via the Selberg integral Eq. (|2^) and calculating 
in this way the fluctuations around the chosen saddle points we find the asymptotic expression for 
the negative moments of the second type: 



= (27r)" 



det {fiilN — H) det (/i2 1 n 
1 



[/ii - ^1*2] 



H) 



(40) 



-nN{l+tJ.^ /2)+iNmrLjp{tJ.) 



n)l 
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which enters the numerator of Eq.(^. Dividing this expression by that presented in Eg. (p6|) and 
taking into account: 



n 



N-n ' 



iV" as TV ^ oo 



we arrive at the announced formula Eq. . 

4 Replica limit. Chiral GUE models 

Let us now briefly consider impHcations of the derived negative moments representations for per- 
forming the rephca hmit n ^ 0. We recall the main steps of the scheme for the positive moments 
("fermionic replica") as suggested by Kamenev and Mezard||ll[, see also Yurkevich and Lerner|]l^. 

Given the expression Eq.(p5|), one takes into account two types of stationary points: the "max- 
imally symmetric" one qi i — q2.i — q'^ as well as all possible sets where exactly one of qi^i and 
exactly one of 52.1 are taken to be equal to q~ , the rest 2n — 2 being equal to g"*" as before. It 
was demonstrated that taking the factors arising from multiplicity of the saddle-points and gaus- 
sian fluctuations around them (Sclbcrg integrals) into accoimt only those two possibilities produce 
leading order contributions nonvanishing in the replica limit. 

Let us stress clearly the bizarre nature of this prescription, as compared with the well-defined 
stationary point procedure for the integer values n = 1,2, First of all, for n = one takes 
only one of two "mostly symmetric" saddle points discarding its partner qi^i = ^2,1 = 9 • At the 
same time, for the positive integer moments the latter produced exactly the same contribution at 
vanishing imaginary part Imfi — > 0. The step is dictated by "causality arguments" pH , i.e. by 
necessity to break analyticity inherent in the positive moments, see [ p9[ . 

Second, the saddle-point sets containing admixture of two q~ contribute now to the same leading 
order as the "fully symmetric" one, whereas for any positive integer n the two contributions were 
different by the factor 1/A^. All this is to remind the reader that presently the replica trick is more 
a kind of art rather than science (or, rather a kind of alchemist's wisdom than regular chemistry). 
For the present author it is however in no way an intimidating characteristics but rather a challenge 
to imagination. 

Let us now turn our attention to the negative integer moments as described by Eq. (^ij) and com- 



pare them to Eq.(35). A little inspection shows that all the factors that make those two expressions 



T-riV-l I 

different are immaterial in the replica limit. For example, lim„^o IlAr-n J- = lini?i^o -r^-iv-n-i = 1 

111 

and the same is valid for [det q\" and other factors. Thus, in the limit n — > the two expressions are 
indistinguishible on the level of saddle point sets and expansions around them. The only essential 
deviation which seems to persist is the difference in the domain of integration, which is half the real 
axis for all the negative integer moments. The latter feature is really dictated by analyticity (or 
causality), which, being a meaningful notion for all negative (but not for the positive!) moments, 
dictates only one saddle point to be operative - that with all q^. 

At the same time, there is no obvious reason why other stationary points should be excluded 
from a consideration in the replica limit. All the experience of working with the replicated ex- 
pressions suggests that saddle-points irrelevant for integer n could be most relevant for n = 0, and 
vice versa. As a distinguished example one can invoke the famous Sherrington-Kirkpatrick model 
of spin glasses where the saddle-points dominating in the replica limit are, in fact, local maxima 
rather than minima of the corresponding functionals. Moreover, formally dominant contributions 
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in that case seem to come from the boundaries of the integration domains, but are discarded as 
"unphysical" in favour of the mentioned maxima, see e.g. discussion in p. 869 of the reference [^^ . 

We therefore suggest that a sensible recipe to perform the rephca hmit for negative moments 
of the characteristic polynomials is as follows: (i) Find an integral representation for the moments 
with help of the Ingham-Siegel-like integrals (ii) Evaluate the resulting integral as a sum over 
the stationary points, starting with the most symmetric set as dictated by analytical structure, 
and adding to it those discovered by Kamenev and Mezard irrespective of the constraints on the 
integration domain. 

To illustrate that such suggestion makes sense beyond the GUE model let us briefly consider 
one more example. This is the so-called chiral GUE introduced to provide a background for 
calculating the universal part of the microscopic level density for the QCD Dirac operator, see 
p3| and references therein. The quantity to be calculated are negative moments of the spectral 
determinant: 



-jvTrJtj 



det 



iJ^ 



iJ 
mlN 



(41) 



where m > is a parameter proportional to the quark mass and, in the simplest case of zero 
topological charge, J is a complex random N x N gaussian matrix, with being its conjugate 
and An being the normalisation constant. 

We relegate the details of consideration of this interesting and important model, as well as its 
close relative - that of non-Hermitian random matrices - to a separate publication p3| and present 
here only a brief account. 

Application of our method based on the use of the Ingham-Siegel type integral Eqs. (0|5|) as 
an alternative to the conventional Hubbard-Stratonovich transformation results in the following 
simple formula: 



dQe 



inTtq 



>o 



det 



(in + niQ 



-N 



detQ 



(42) 



which is exact for arbitrary N > n. Here Q is a positive definite n x n Hermitian and the constant 
is given by: 

1 



In the thermodynamic chiral limit one considers m ^ , N oo but keeping the product mN = ^ 
fixed. This results in reducing the above expression to: 



.N 



dQe-fTr(Q+Q-^) 



Q>0 



detQ 



(43) 



This should be compared with the corresponding formula for the positive moments p2 

I'^iim) = A^/dJdjte-^Tr,;ti^^V™i^ ^J 
J \ iJ^ mlN 



(44) 



oc 



dfi{U) exp{--Tr (u + U-^)} 



Uin) 



where the integration goes over the unitary group U{n) 
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The correspondence between the two integrals is very similar to the GUE case discussed by us 
earlier in this paper. Considering x — > oo as the parameter justifying the saddle-point approxima- 
tion one finds the saddle-point sets of both integrands coincide: they are given by matrices with the 
eigenvalues ±1. Again, in view of the constraint Q > any negative integer moment is dominated 
by the most symmetric set with all eigenvalues being equal unity, whereas positive moments are 
given by the sum of contributions of many such sets. The Kamenev-Mezard-Yurkevich-Lerner lim- 
iting procedure which uses the most symmetric -1-1 configuration as the reference point was shown 
to produce sensible results for the integral Eq.(Q) Taking into account that the difference 

between the two integrands is immaterial in the replica limit, we again arrive at the conclusion 
that we can not help but adopt the same scheme for proceeding from the negative moments. 



5 Implications for the supersymmetry method 

A curious point which is worth mentioning is that for the case of positive moments the use of the 
Hubbard-Stratonovich transformation is, in contrast, very natural and effective. Attempting to use 
our method for that problem encounters with the difficulty of dealing with diverging integrals. The 
latter hide, in essense, necessity to work with higher derivatives of the S— distributions. All this 
is suggestive of a certain duality between the two methods: working with Grassmann integrations 
requires the use of the Hubbard-Stratonovich transformation but dealing with commuting vari- 
ables is much facilitated by avoiding that route. This observation has certain implications for the 
supersymmetry method treating both types of the variables on equal basis. From this point it is 
of independent interest to try to treat those two differently, employing the Ingham-Siegel integrals 
of second type, Eq.(^5|), for the commuting ones. Such a treatment for the simplest nontrivial 
example of the generating function containing simultaneously both positive and negative moments 
of the characteristic polynomials of GUE matrices is presented below . We will see that our method 
indeed works and has certain advantages in comparison with the standard route. For example, the 
so-called "boundary terms" which usually arise due to a singular nature of the transformations 
along the standard route would not appear in our method. 

We are interested in calculating the following generating function: 

CN{HB,f^F) - ( r, / , . . X ) (45) 

\ ZN(p-lb)^N(P2b) I GUE 

where fiB = diag(/iif,, /ijb) i /*f — diag(/^i/, //2/) and Im(/iib, fj,2b) > 0. That generating function is 
obviously an analytic one in {^if,^2f) plane. It turns out to be technically convenient to change: 
fiif —i^iif , fJ.2f ~~i^^2f when performing the ensemble averaging, and restore the original 
generating function by a simple analytical continuation. 

To calculate the average we first use the standard " supersymmetrisation" procedure and rep- 
resent each of the two characteristic polynomials in the denominantor as the Gaussian integrals, 
Eq.(^, whereas those two in the numerator are represented as the Gaussian integrals over the 
anticommuting (Grassmannian) A^— component vectors Xi,Xi^ and X2,X2 j ^-S- lUl- The 
ensemble average is then easy to perform and after a straightforward manipulations the generating 
function can be brought to the following form: 

CN{fib,fif) oc J dxidxl J dx2dxlexpl^^(^nifXiXi+ fJ'2fX2X2) + -^Tt(q%^^ 

d^Si yd2s2exp|^/xiStSi-^/x;stS2-^Tr(gsiQi3^)| (46) 
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where we introduced the 2x2 matrices 

n ( x\xi x\x2 

V X2X1 X2X2 



Si (g) si - S2 < 





Qi 



s|Si s|s2 

S2S1 S2S2 



X2 



and L = diag(l, — 1). 

Now, following the main idea outlined in the beginning of the present section we employ the 
Hubbard-Stratonovich identity: 



exp 



1 

8iV 



— TrQl 



cx / dQ F exp 



1 



2^JN 



X2 



(47) 



where the integration goes over the manifold of 2 x 2 Hermitian matrices Qf = 



qii qi2 

9l2 122 

Changing the order of integrations and exploiting the above identity one performs the (Gaussian) 
Grassmannian integral explicitly in a simple way and brings the expression to the form: 



(48) 



dQpe 



-^TrQl 



(i^Sifi^S2exp<j ^/xiSjSi - ^^28^82 - ^Tr [QbLQbL 



det 



Ml/ " 7^911 



1n-B 



J^qi2lN 



where we introduced the N x N matrix B 



1 

2N 



M2/ 



Si 



1 

/n 



922 



B 



Observing that TrB" oc Ti{QbL)" for any integer n, and introducing Ai/,A2/ as two (real) 
eigenvalues of the (Hermitian) matrix Qp ^ jlp — '^^^ '^^^ satisfy oneself that the determi- 
nant factor in the integrand of Eq.(p^) is equal to: 



(detg^)^- 



■det 



det 



A2/I2 - -^QbL 



Next step is to evaluate the integral over 81,82. Since the integrand depends on those variables 
only via the matrix Qb the calculation can be done exactly the same way as elsewhere in the 
paper. That involves the (5-function representations and the subsequent Ingham-Siegel second type 
integration, see Eqs. (f3|,E5 16). The resulting expression is given by: 



CAr(A&,A/) oc 



dQp (detQ^) 



^,7V-2^_iTrQ^, / dg^dctg^-^exp 



1 /~"\^ i r 

Tr [QbL j + -Tr \jxbQe 



8N 



X det 



det 



^^2/12 - -^QbL 



16 



Further steps involve: (i) changing: Qf vNQp, Qb '^NQg (ii) diagonahzing QbL = 
Tdiag{pi,p2)T~^ , where pi > 0,P2 < (iii) integrating over the "hyperbohc" manifold of the 
pseudounitary matrices T in the same way as it is done in Eqs. ( ]29| , ^0[ ) (iv) shifting the integration 
over Qf to that over Q^' and introducing the (real) eigenvalues qi, q2 of the latter matrix (and the 
corresponding unitary matrices of eigenvectors Up) as new integration variables, (v) performing the 
integration over the unitary group Up by the standard Itzykson-Zuber-Harish-Chandra formula, 
see Appendix C. Continuing finally analytically as (/ii/,/i2/) ^ (*Aii/i*/^2/) we arrive at |^ 



N / 2 I 2 \ 

CN(fib,tlf) cx e-(^i/+^^/) 
dqi 
dpi 



{qi - 92) 



dq2 

O 

dp2 



{Pl +P2) 
{VIV2Y 



exp{~iV [Lif{qi) 
exp{-A^[£if,(pi) - 



- £2/(92)]} 

>C2fc(p2)]} 



(50) 



where 



£1/(9) = 



2" 



2^ 



- In? 
i/iihp- Inp 



C2s{q) = \q' 



«M2/g - 



In q 
lap 



(51) 
(52) 



So far all the formulae were exact for any integer N > 2. Let us stress an apparent symmetry 
of the resulting equation with respect to q and p variables, apart from the anticipated difference 
in the integration domain. Consider now the limit N 00, where the integrals are dominated by 
the saddle-points: 



Pl 



t4b 



qif 



P2 



1 r 
2 



i/-i26 + 



(53) 



M2fc 



Here we took into account the restrictions Repi 2 > and also set Im/i(2.2)6 — 0. 

In fact, in the limit N 00 we are interested in the "local regime" : M(i,2)b = /^i ^ i M(i,2)/ 
/i ± ^ where < 2 and uji, , ujf — 0{1/N). We find it convenient to define: 



M(l,2)/ 



— 2 sin ( 



(1.2)/ 



Ai(1.2)h 



= 2 sin ( 



(l,2)h 



Then Ubj = Mi(h,/) - ^^2{b.f) ^ i<t>i{bj-) ~ '/'2(fcj))2 cos 0, where 2 cost 
proportional to the mean spectral density. Thus we have: 



= a/4 - ^2 = 7rp(^) is 



916 



P2b = 



7rp(/i) 



91/ 



52/ 



^/ 

7rp(/i) 



The saddle-point values, Eqs.(|5^) are given by: 

P2 



Pl = e"^" 



„— 2</>2b 



qf = e^^^'^'i^ 



■^One should take into account a symmetry of the integrand with respect to the transformation: qi 
simphfy the expression. 
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which shows that: 



9i 



- qt = e^">''' (l - e±'('^i/-'^2/)) ^ 0{\/N) 



it 



- = ±2 C0S(/.(i,2)/ + 0{1/N) « ±^^pi^^) 



The latter relation makes it clear that only two saddle point sets {q^ i q^) give dominant contribu- 
tion whereas the other two: {q^ ^ q^) are suppressed due to the factor {qi — 92) in the integrand of 
Eq.®. 

Expanding around the relevant saddle points and picking up the leading order contributions 
we obtain after a standard set of manipulations: 



JN{<pii-,p2b) . 



ifJ-lb - /i2b)(/il/ - M2/) 



„iAf(0i/-02/) (pi'Plf _ pi't'lb) 



-i<t>21 



-i<p2b\ 



CX 



UJbUJf 



_g-lJV(0iy-</,2/) 

1 



(54) 



The latter expression coincides with the particular case of the result obtained in p3| by a different 
method. 

The method based on the combination of the Ingham-Siegel and Hubbard-Stratonovich trans- 
formations proves also to be a direct way to derive the mean density of eigenvalues for the chiral 
GUE ensemble discussed in the previous section. Let us outline the corresponding calculation for 
the (quenched) case of zero topological charge. 

The generating function to be calculated is given by (cf. Eq.(B^)): 



JN{fnf,mb)(x / dJdJ^e 



-nTtPj 



det 



det 



i,P rublN 



(55) 



where mf,mb > 0. 

Representing the determinant in the numerator and that in the denominator as the Gaussian 
integrals over anticommuting and commuting variables respectively, we easily perform the ensem- 
ble averaging, and then employ the Hubbard-Stratonovich transformation for the terms quartic in 
Grassmannians whereas using the Fourier representation of (5-function for their commuting coun- 
terparts. This allows to integrate out Grassmann variables exactly and after the Ingham-Siegel 
integration to arrive to the following expression: 



l!S[{mf,mb) oc J dqdq*{q*q)^ ^exp— ^ 

f°° f°° N ( 

X J dpi J dp2(piP2)^"M9*9-PiP2)exp-— |pip2 



(56) 



TOfc 



(Pl +P2)| 



Ghanging now to the polar coordinates R, 9 in the complex g, q* plane and next to the variable 
t — P1P2 one can perform the 6 and pi-integrations explicitly finding: 



lN{xf,Xb)(xe " 



dr 



dt{rtf-^ (r - t) e-^(''+*)/o(2x/ri/2)i^o(2a;faii/2) (57) 
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where x^j = Nnibj/A and /o(2),ifo(^) stand for the modified Bessel and Macdonald functions, 
respectively. In the most interesting chiral hmit we set N —f oo while keeping x^j fixed. The 
saddle-point values are obviously given by ?' = t = 1 and expanding around them up to the first 
non- vanishing term finally yields: 

Ichtrai{xf,Xb) cx [x f Ii{2x f)Ko{2xb) + XbIo{2x f)Ki{2xb)] (58) 

in full agreement with known results, see [^3[ 



6 Conclusions and Perspectives 

In the present paper we suggested a systematic way of evaluating negative integer moments of 
the (regularized) characteristic polynomials. Using the standard representation of those moments 
in terms of the Gaussian integrals as the starting point we found a route avoiding the use of the 
ubiquitous Hubbard-Stratonovich transformation. Instead, we advocated the exploitation of the 
matrix integral Eq.(^5|) similar to that considered long ago by Ingham and Siegel, Eq.(0). The 
advantage of the procedure is that the emerging structures are attractively simple and, in essense, 
very close to those derived earlier for the positive moments. We evaluated the resulting integrals 
in the limit ^ oo by the stationary phase method and extracted the leading asymptotics. The 
limiting value of the correlation functions for the negative moments is presented in Eq.(Q) and 
expressed in terms suggesting universality. In fact, by using the character expansion method we 
were able to reproduce the formula (jj) for unitary random matrices|Q. We therefore conjecture 
that it should be equally applicable, mutatis mutandis, to the behaviour of the Riemann zeta 
function in the close vicinity of the critical line. 

Our analysis may raise a few questions which deserve further discussion. First of all, one may 
wish to know if it is possible to arrive to the same representations via the standard (Hubbard- 
Stratonovich) method. The answer is of course affirmative as is demonstrated in the Appendix D 
on the simplest non-trivial example. The general case can be treated along the same lines. We 
however insist that such a way is hardly natural for the present problem and, in fact, obscures the 
simple structures arising. 

This fact has important implications for the " supersymmetry method" of treating generating 
functions involving simultaneously characteristic polynomials both in positive and negative powers. 
We presented a few simple non-trivial examples of such calculations in the text. In fact, we found 
it natural to treat "bosonic" and "fermionic" auxilliary integrations differently, which was a serious 
departure from the standard spirit of " supersymmetry" . 

One should mention that the correlation functions calculated above can be also calculated in the 
framework of the standard Efetov's supermatrix approach. However, present method, to our mind, 
has some technical advantages. For GUE calculations, for example, the two terms appearing in 



the final result (54) on equal footing have quite a different origin in the Efetov's method. Namely, 
one of them appears in the form of the so-called " anomalous" , or " boundary term" . When trying 
to calculate the correlation functions of the higher orders the boundary terms turn out to be quite 
difficult to keep track of, making the cr— model calculation a daunting job. 

In fact, the authors of announced the result for the most general expression of the ratios 
of integer powers of characteristic polynomials . They arrived to it by a method generalizing that 
suggested in the paper by Guhr[^ and employing an analog of the Itzykson-Zuber integral over 
a supermanifold. Our method provides an alternative way of derivation of the higher correlation 
function, the corresponding calculation will be presented in a detailed form elsewhere [Eil 
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For the case of chiral GUE ensemble the apphcation of the standard Efetov's method for 
the quenched model requires quite a cumbersome calculation [|5|, whereas along our route the 
calculations are rather short and elementary. Moreover, we found that the same method 
proved to be applicable for the most interesting " unquenched" situation when the Efetov's method 
encountered with unsurmountable difficulties. 

Finally, it is interesting to explore if the Ingham-Siegel integrals and their natural generalisa- 
tions could provide a serious alternative to the Hubbard-Stratonovich transformation in the whole 
class of problems in the domain of random matrices and disordered systems. To this end two 
aspects are worth mentioning: (i) the Ingham-Siegel integrals of both types are known for all sym- 
metry classes, see Appendix A and (ii) performing the saddle-point calculation directly on the level 
of Eq.(p8|) yields the standard non-linear cr— model representation for the negative moments. 
Further work along these lines is under the wayp3, 44 but a general affirmative or negative answer 
to the questions requires more efforts. 
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APPENDICES 

A Calculation of the integral Eq.(15) 

Our goal is to calculate the integral 

In,N{Qn) = J dF„e5Tr(^^.Q.) [det - ^lln)Y^ (59) 

where both and Qn are Hermitian n x n matrices. First notice that the integrand is invariant 
with respect to the unitary rotations F UFU~^, hence the result of the integration can depend 
only on the eigenvalues of Q. Then, following |3^, ^ one can take Q to be diagonal from the very 
beginning and separate the first eigenvalue from the rest: 

Q = diag(gi,g2, ••■,9«) = diag(9i, (3„_i) 

Accordingly decompose the matrix Fn as 

Fn^(^ "^J' ^ , dF„ = dfndtUtdP^., (60) 

where = (/|]^, f^i, f*i) is a n — 1 component complex vector. 
Next step is to use the well-known property of the determinants: 



det 



(Pn ~ = det (^K-i - Mln-i) (^fii F„_i - /il„_i 
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which gives: 

l7i,N{Qn) 



J J — OO 



dct ( Fn-l - filn-1 
1 



(61) 



/ii - M - ft 



TV 



The last integral over /n is evaluated by the residue theorem taking into account Re/i > 0, the 
result of the integration being: 



2TTi 



%i)(^) exp^^gJ^ + ft 



N-l 



Fn-l — Ml"- 



(62) 



where 9{x) — 1 for x > and zero otherwise and we assumed > 1. Now the gaussian integration 
over df^df can be easily performed, yielding the factor: 



so that we arrive at the recurrence relation: 



In,N{Qn) = ^ ^^^Zi ( ) ^^(gi)e*^''^/«-l,JV-l(Q„-l) 



T{N) \2 

which immediately produces the desired formula 



J-n,NW7i) — ' T^Ty 



-n^fe)det 



uTtq 



(63) 



(64) 



assuming N > n. 

In fact, the derivation is straightforwardly repeated for the case of real symmetric matrices F„ 
and Qn- The recurrence relation in that case is: 



■n ^ , 



T{N) \2 



o9i 



N- 



which yields the result: 



I-%{Qn) = r '^l n 0(<z,)det 



LiTrQ 



(65) 



(66) 



for N > 



B Properties of the matrices Ql 

In this Appendix we consider the manifold of 2n x 2n matrices Ql — QL where Q^ — Q > and 
L = diag(l„, — 1„). In fact, this set of matrices is closely related to the object known as a regular 
matrix pensil |^6|. 
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We begin with proving that all eigenvalues of such non-Hermitain matrices are real and half 
of them positive, the rest being negative. In doing this we can safely assume that all eigenvalues 
are different since matrices with degenerate eigenvalues form a manifold of lower dimension and as 
such will not contribute when we integrate over the whole manifold of Ql- 

The characteristic polynomial for the eigenvalues q of the matrix QL can be written as: 

det (gl2„ - gi) = detgi/2 (gl2„ - Q^^^LQ^^^^ Q^^/^ ^ (^^^^^ _ q1/21q1/2^ 

where we used that Q^/^ > is a nonsingular Hermitian matrix. Then all eigenvalues of Qi^ 
coincide with those of the Hermitian Q^I^LQ^I'^ and therefore are all real. Moreover, the number of 
positive and negative eigenvalues of any Hermitian matrix H stays invariant under transformations 
ii — !■ T^HT, where T is an arbitary nonsingular matrix ||4q| . We arrive at the conclusion that the 
number of positive and negative eigenvalues of Q^^^LQ^^is the same as that for L, thus proving 
the statement. 

Let qj be an eigenvalue of Ql and denote the corresponding (right) eigenvectors as e^: 

(qij ej = qjBj , e] {LQ^ = q^e] 

Multiplying the first of these relations with e\L from the left and the second relation with Lej 
from the right we have: 

el {LQLj Gj = qjelLej ^ qkelLej 

showing that e^Lej = SjkqjGjLej 

Now, L is a unitary matrix, hence LQL > so that qje'jLej > 0. Introduce now the "normal- 
ized" eigenvectors = ej/^ sgn{qj){ejLej), where sgn(2;) stands for the sign function. Then it is 
easy to see that 

e'jLej = sgn(qj) and e\LQLej = \qj\5kj 

Further introduce the matrix T whose columns are vectors e-,- for j = 1,2,..., 2n, and consider 
Tl = L. It is immediately clear that 

TlT = diag(sgn(gi), ...,sgn(92n)) and TlOLT = diag(|gi|, |g2ri|) = TLrdiag(gi, 92™) 

The vectors are obviously linearly independent, hence the matrices and T are nonsingular. 
This immediately shows that the matrices QL can be diagonalized by a similarity transformation: 

QL = fdiag(gi,...,g„)f-i 

where T satisfies: 

f^LT = diag (sgn(gi), sgn(92«)) 

The last matrix is essentially L up to a permutation of its entries on the main diagonal. This 
completes the proof. 



22 



C Evaluation of the integral Eq.(|5D|) 

To evaluate the quoted integral one needs to employ an explicit parametrisation of the matrices 



To e 



U(n,n) 



^ U{n)xU(n) 



. We follow the paper pa] where it was suggested that the following parametrisation 



is especially convenient: 



i 



\/T+W 



in terms of complex n x n matrices t , . The reason for such a choice is dictated by an especially 
simple form of the integration measure: 

Next step is to diagonalise F with help of two unitary rotations: ft — u'^tuji where ua,r G 
U{n) and f = diag(Ti, ...,r„), so that i = u'^t'^ua. It is convenient to write the modulus and 
the phase of Tk exphcitly: Tk = sinh V'fcc"^'', <7/'<oo,O<0<27r,A: = l,...,n. 

The matrices Tq take the form: 

cosh ip e"^ sinh -0 \ / 



( 1" 


)( 







and Tq"^ is correspondingly given by 

--1 



sinh -0 cosh ip 



UR 



^0 — 



u, 



u 





-1 

R 



UA 

iiR 



cosh ip — e*''' sinh %l> 
e~*'^sinh-0 cosh^/; 

One can straightforwardly calculate the integration measure in the new variables and find: 

n n n 

2 

fc=l ki<k2 k=l 



didP oc ]^sinh2V'fe (sinh^ -0^^ — sinh^ i/j^^)^ ]^ rfV'fe c^^fe x (67) 

fcl <fe2 



n 



n 



« n ('^*=i ^ '^''^2)^ n c^'^fe d(pk X dfi{uA) dfi{uR) 

ki<Ck2 k—1 



(68) 



where dfi{uA,R) are normalised invariant measures on C/(n) and we introduced: = cosh2'!/'fe G 
[1, 00) as new variables. 

Now we can use the cyclic permutation of the matrices under the trace sign to rewrite the 
expression in the exponent of Eq.(pO) in terms of the introduced variables as follows: 



Tr 

Tr 
Tr 



Pi 



P2 



f-i 



cosh ip e**^ sinh ip 
"*'^sinh?/; cosh?/; 



Pa 



Pe 



cosh Ip — e*"^ sinh tp 
— e^*"^ sinhf/) cosh?/; 



^ f 2^ ^ f 2^^ 2^^ 

Pa I Ml cosh Ip ~ ^2 sinh ?/; 1 + Tr ( ^3 cosh ip — ^li sinh -0 



= ^(a^i + Ai2)Tr (Pa + Tb) + i(Aii - M2) (TrT4 cosh2?A - TrP^ cosh 2^ 

where we introduced matrices Pa — uaPiuJ^^ , Pb — UBP2'Ug^ having the same eigenvalues qi — 
diag(qi^i, ...,9i,„) and q2 = diag(q24, ■•■,92,n) as the matrices A,2- 
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We see that the integral of interest is expressed now as: 



cx / TTdAfc A2(A)e5^('^i+''^)Tr(«i+92) (jq) 
•^1 k 

X / d^(uA)e5^(''i-f2)Tr«APifl^^A /" d^(y^)e-5^(A'i-M2*)Tru«P2«^U 



where we used A = diag(Ai, A„) and the symbol A(A) for the corresponding Vandermonde 
determinant. 

Two integrals over the (normalized) Haar measure on the unitary group are given by Harish 
Chandra-Itzykson-Zuber formula[|^ : 

where in our case (3 — ^N{fii — /I2) for the first integral, and for the second one /3 — /?. This 
gives: 

/(M,A,P2) oc \ f^(p.+M;)Trte+,-.) 

A(gi)A(g2) 

r r n ^ r n 

X / TT^^Afcdet e*^(''i-^2')^'=''i'' det e^^^^^i-^^'^^'-'^.' (71) 
J I I Jfe,i=i L Jfc,;=i 

The last integral can be easily calculated by expanding each of the two determinants as: 

N 

det = ^(-l)[^lexp{±/3^ A.g.J 

[S] k=l 

where [S] stands for a permutation ri,r2, ...,r„ of the index set 1,2, ...,n. Product of two such 
expansions can be integrated term by term, and the integrals are convergent due to Re/3 > , qi^i > 
, q2.i < 0. This gives: 

/°°n^A.det K^--]:,^,det [e-'^^--];,^, = [-/?]-" ^ (-1)'^^'-^'^-' 

where [Si], [82] are two independent permutations (ri, r„) and (^i, In) of the index set 1, 2, n. 

Clearly, one can restrict the summation to be taken over the relative permutations of the two 
index sets and multiply the result by n\. The exponential above is invariant with respect to any 
index permutation, so it can be taken out of the summation sign and one recognizes the so-called 
Cauchy determinant: 

1 V A(gi)A(g2) 



in the remaining sum. Collecting all the relevant factors together we arrive at the final formula: 



I{M,Pl,p2)^[-^{pl-^l;)Y 



up to a constant factor, which can be fixed by normalisation in the corresponding equations. 
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D Negative moments by the Hubbard- Stratonovich trans- 
formation 



Let us satisfy ourself that the standard Hubbard-Stratonovich transformation over the hyperboHc 
manifold p8|, pl| produces the same formula Eq.([3l|). We concentrate on the simplest nontrivial 
case n = 1 for the sake of clarity. That case was used for a pedagogical introduction into the 
Hubbard-Stratonovich method in the author's lectures in the book and the notations mainly 
follow those lectures. 

Our starting point is Eq.(|27|) for n = 1. We introduce the matrices A ~ L^^^QL^/^ so that 
Tri2 ^ TtQLQL and TtA = S\Si - SjSa , TrAL = S|Si + 8^83. Remembering /^i,2 = M ± 



{lu/2 + iS) we can express all terms appearing in the exponent of Eq.(27) in terms of A: 



[Si - ^^28^82 ) = '-^iTl■A +1(6- tiu/2)TvAL 



The Hubbard-Stratonovich transformation is the identity: 



— Tri^ --(S- iuj /2m AL 
8N 2^ ' ' 



oc e 



exp 

Ar(i5-iw/2) 



(73) 



dQ exp 



TrQ2 - iN{& - iw/2)Tr((5L) - -TrQi 



Despite looking as an innocent gaussian integral the identity is very nontrivial, since the con- 
vergency arguments force one to choose the following "hyperbolic manifold" of the matrices Q as 
the integration domain: 



Q 



T- 



T 



T ■ 



cosh 6* e^-^sinhe* 
e-*'^sinh6' cosh 6* 



(74) 



Pi 

P2 

dQ oc [pi — P2f' sinh 29 dpidp2 dO d(f>, 

-00 < Re(pi,2) < 00 , Im(-pi,p2) >0,0<6'<oo,0<(?!)<27r 

A detailed discussion of the convergency problems and of the above identity Eq.(|7^) can be 
found, e.g., in the book and in the mentioned lectures 

Substituting such an identity back to Eq.(27) and changing the order of integrations over Q 
and Si 2 we see that it can be processed as follows: 



^N(S~ioj/2f 



dQ exp <j -y TrQ^ - iN{5 ~ iLo/2)Tv{QL) 



(75) 



j rfSrfSl^expj; 



X / dSrfStexp<{ ^^StLS- ^StLi/2QT^Li/2s 



^N{S-iuj/2y 



-N 



dQ det(^/„-Q) exp<j -yTrQ^ -i7V((5-itj/2)Tr(QL) 



dpi / dp2- 



[Pl - P2f 



[(Ai-Pi)(M-P2)] 
d9 sinh 26* exp | ~iN{S - iuj/2)Ti 



exp 



N , 



-{pI+pD 



PI 

P2 
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where we introduced the notation S = (81,83) and used TiQA = S'' L^/^Q'^ L^/^S. 
Using the exphcit parametrisation for T it is easy to verify that: 



Tr 



P2 



(pi — P2) cosh 20 



which allows one to perform the integration over 9. The expression above is therefore reduced to 



dpi / dp2 



P\ exp j-f + pI) - ^NiS - zc./2)bi - p^) } (76) 

[{fl - Pl){H - P2)\ I ^ 



S-iuj/2 

Remembering Impi < , lmp2 > we can use the identity 



[(m-Pi)(m-P2)] ^ oc 



dqi 



dq2{qiq2) exp{i [(^ - pi)qi - {fi - ^2)92]} 



and further rewrite the integral Eq.(^) as: 

OC — ■ — — I dqi J dq2{qiq2)" 



^-lp«M(9i-92) 



(77) 



N 



{Pi+pD - iN{5 -iuj/2){pi -ps) - i{piqi -^292) 



5 — iu}/2 Jq 
X J J d.P2{pi-P2)exp '^ ^ 

The next step is to use the chain of identities: 

j dpi y"dp2(pi-P2)exp|-y(p^+p2)-«A^('5-«w/2)(pi-p2)-i(pi'?i-P2'72)| (7 



9 9 \ f 1 

' ^ + ^)^'^P|"^(«i+9|)-('^-»w/2)(gi + g2)-iV((5-iw/2)^ 
and observe that 



0^(91-92) 



g^jn<Zl,.2)0c(^— + — 



pV('Jl-'j2) 



The last formula allows one to integrate by parts over qi , 52 and in this way to get rid of the 
derivatives. The boundary terms vanish for > 1, the application of the operator 

+ to 

(9192)^"^ produces the term (<Zi'72)^^^('?i + 92) and the resulting expression concides with that 
given in eq.(pl|) for n = 1. 

For n > 1 the equivalence can be shown along essentially the same lines, but calculations 
become cumbersome and require the use of the procedure similar to that outlined in the Appendix 
C. 
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